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Abstract 

The z-zeros of the modified Bessel function of the third kind K v {z), also known 
as modified Hankel function or Macdonald function, are considered for arbitrary 
complex values of the order v. Approximate expressions for the zeros, applicable in 
the cases of very small or very large \u\, are given. The behaviour of the zeros for 
varying \v\ or argi^, obtained numerically, is illustrated by means of some graphics. 
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1 Introduction. 



The relevance of Bessel functions to the solution of a great variety of problems 
in Physics and Engineering is widely known and does not need to be stressed. 
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In particular, zeros of the different kinds of those functions are closely related 
to energies of bound or resonant physical systems. For instance, the i/-zeros of 
the Hankel function H u (z) for positive values of z determine the poles of the 
amplitude of scattering of various kinds of waves by spheres and cylinders |2Uj . 
Almost four decades ago we considered [10] the zeros of the modified Bessel 
function of the third kind, K u (z), also known as modified Hankel function or 
Macdonald function. Information about the z-zeros of K v (z) for positive order 
v was found in the book by Watson [33]. However, our interest, motivated by 
a quantum mechanical problem, was on z- zeros when v assumes purely imag- 
inary values. For that case it was already known [T2"]n%] that K v {z) presents 
an infinity of zeros, all of them located on the positive real semiaxis of the 
complex z-plane. We showed that these positive zeros, labelled according to 
decreasing values, form a sequence tending to be a geometric progression with 
ratio exp(— 7r/|z/|). We gave, also, asymptotic approximations to the position 
of the largest zeros for large values of \v\. Further discussions of the properties 
of the zeros of K u (z), v pure imaginary, have been carried out by Laforgia [21] 
and by Dunster [7]. Zeros of K v (z) with complex v have been less studied. In 
view of their application in the analysis of electromagnetic wave propagation 
in anisotropic inhomogeneous conducting media, Nalesso [20] has considered 
the i/-zeros of K v {tv) for fixed positive values of t. But, to our knowledge, a 
discussion of the z-zeros of K u (z) for arbitrary complex v is still lacking. 

In recent years, stimulated by review papers by Lozier and Olver [2"3"f2~S] and 
by the DLMF project, a considerable number of new algorithms for the com- 
putation of special functions have been published. Concerning the Macdonald 
function, Refs. [9lH3fl4lll5lfl 6 17 ,29] are to be noticed. As a general rule, algo- 
rithms lose accuracy in the vicinity of the zeros. Hence the interest of having a 
knowledge, at least approximate, of the location of those zeros. On the other 
hand, K u (z) presents only a finite number of z-zeros in the Riemann sheet 
— 7r < argz < it if v is real, whereas an infinity of zeros occur if v is pure 
imaginary. It has therefore seemed to us interesting to study the evolution of 
the z-zeros of K u (z) as the modulus and/or argument of v are continuously 
changed. 

In what follows, we consider only values of v in the quadrant 
< argz/ < tt/2, 

and z is assumed to lie in the principal Riemann sheet, — n < argz < tt, for 
large values of v. The symmetry relations 

K- V (z) = K„(z), K v {z) = Kjz), (1) 
allow to extend our results to values of v in all quadrants. On the other hand, 
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the relation [Q Eq. 9.6.4] 

K v (z) = Ime^H^ (ze?™) , (2) 

valid when K and H^ 1 ' are analytically continued to every Riemann sheet, 
proves that the pattern of z-zeros of the Hankel function HjP(z) results from 
that of K v {z) by rotation by an angle of +7r/2 around the origin. 

We discuss in Section 2 the behaviour of the zeros for decreasing values of \u\. 
The case of large values of \v\ is considered in Section 3, where asymptotic 
approximations to the zeros are given. Numerical results for the zeros of K u (z) 
with moderate \v\ are presented graphically in Section 4. Finally, the particular 
case of real v is commented in Section 5. 



2 Zeros of KJz) for \v\ < 1. 



The particular case of pure imaginary v revealed [10] that, as \v\ goes to zero, 
all zeros of K v (z) approach the origin. Therefore we start by assuming \z\ <C 1 
in the expression of K u (z) as a sum of two convergent ascending series P, Eqs. 
9.6.2 and 9.6.10] 



sin UTT 



( {z/2)- v ~ (z 2 /4) k {z/2Y ~ (z 2 /4) k \ 
\T{l-v) to k\{l-v) k r(l+i/) to M(l+v)k) ' 



valid for v different from an integer. Obviously, the zeros of K v (z) should 
satisfy the relation 

~ {z 2 /A) k 
'l+v) hk\(l- 



r(l + u) un k\(l-u) k 

(z/2) 2u - y ' — - — (4) 

u ' ■ r(i-i/) g (z 2 /4) k ' 



that, in the case of being \z\ <C 1, can be approximated by 

<^"-^( i +2(iV 2+o <4 N<<L (5) 



By taking logarithms, one obtains, as a first approximation, 
log |*| +iargz ~ + log 2 + — log(r(l+i/)/r(l-i/)) + — ^, (6) 
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where n could in principle take the values n = 0, ±1, ±2, . . ., although, accord- 
ing to that obtained in the case of v being pure imaginary, actually existing 
zeros correspond to 



n = 1,2,3, .... 



These values of n will be used as a label to characterize each zero in the form 
z n . Ignoring, in the crudest approximation, the last term in the right hand 
side of Eq. (jSJ) , one has 



\og\z n \ ~ M (-(-niTi + ±log(r(l+v)/r(l-v)fj +log2, (7) 
argz n ~3 (-{-mri + \ log(r(l+i/)/r(l-i/)) 

The behavior of the zeros z n as \v\ decreases can be immediately obtained 
from these expressions by using the series expansion [H Eq. 6.1.34] 

1 + 71/ + c 3 u 2 + ... (9) 



r(i+i/) 



to obtain 

r(i+i/) , , 

1 - 2jl> + 2 7 2 z/ 2 + . . . , (10) 



r(i-iz) 

where 7 represents the Euler constant, 7 = 0.5772156649.... By retaining 
terms up to 0(|z/| _1 ), it turns out 



log|z n | ~ -nix — - + log 2 - 7, (11) 
\ v \ 

Mis . , 

arg ^„ ~ -n7T — -. (12) 

r I 

These two equations allow us to obtain a picture of the behaviour of the zeros 
of K v (z) os\v\ tends to zero. For instance, let us assume that \v\ decreases with 
arg v fixed. If arg v = 7r/2, the zeros z n approach the origin along the positive 
real semiaxis, as described in Ref. [10J. If < argz/ < n/2, \z n \ tends to zero 
whereas argz n decreases without limit. In other words, the zeros approach the 
origin spiraling clockwise infinitely. According to Eq. (fl2|) . most of zeros (or 
even all of them, for sufficiently small \u\) occur outside the principal Riemann 
sheet. For arg v = 0, Eqs. (fTTT) and ([121) are no more valid because, as we will 
see later, the possibility exists of z n going to infinity for certain values of v. 
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The behaviour of the zeros in this case is much more complicated and will be 
described in Section 5. 



3 Zeros of K v {z) for \v\ > 1. 



We have already mentioned that Eq. (j2J) allows one to obtain the zeros of 
the Macdonald function from those of the Hankel function, and viceversa. 
Previous research, by other authors [31I11II19|I20[25II27] . on the i/-zeros of the 
Hankel function shows that, for large values of \z\, they occur at v ~ z. 
Asymptotic expansions of H^ 1 ' valid in the transition region should then be 
used to obtain approximate expressions of its £-zeros for \u\ ^> 1. From Eqs. 
9.3.23 and 9.3.24 of Ref. [T], by using Eqs. 9.1.3 and 10.4.9, a convenient 
expansion can be written, namely 



+ 0^1/3) „ * e -/3 Ai(-2 1 / 3 e l27r / 3 ^) £ f k (6)u' 2k / 3 

V \ k=0 

O 1/3 oo \ 

e^/ 3 Ai / (-2 1 / 3 e 12 ^^) V g k (6)v- 2k/3 , (13) 
v 1 k=o J 

where fo(9) = 1, the other functions fk(0) and gk{0) being given in Eqs. 9.3.25 
and 9.3.26 of the same reference. A more powerful asymptotic expansion of 
exists [U Eq. 9.3.37], but its coefficients are much more complicated and, 
although reversion of that expansion could be done by means of a procedure 
due to Fabijonas and Olver [8], we find preferable, for the sake of simplicity 
and transparency, to work with Eq. (fl3l) . A simple inspection of that expansion 
shows that, for large \u\, it vanishes for values of 6 verifying 

9 S = -2- 1 / 3 e^ 3 (a s + e a (u)), s = 1, 2, 3, . . . , (14) 



where a s represents each one of the zeros of the Airy function, given in Ta- 
ble 10.13 of Ref. PQ, and e s (u) = 0(u~ 2 ^ 3 ) is to be determined. Trying an 
expansion of the form 

00 

£» = 5> SJ ^'/ 3 , (15) 



and using Taylor expansions of the Airy function and of its derivative around 
a s , it is straightforward, with the aid of a computer algebra software package, 
to check that the right hand side of Eq. (TIB]) cancels out for 
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^ 3= 12^00^ (479 ^- 40) ' 

asA = 16170000 2_1/3 e ~' M3 a " (2 ° 231 G ' + 551 ° 0) ' 6tC - 
Approximate values of the z-zeros of HjP(z) are then given by 

z s ~ (u + e s u^ 3 ), k|>l, s = l,2,3,.... (16) 



Reversion of this relation between z and z/ would provide an approximation 
for the i/-zeros of H)p(z) when |z| is large, 



v s ~ 2 + 2- 1 /V i27r / 3 a s z 1 / 3 + — 2 1 /V 2 ' r / 3 a= z^ 3 - —(a 3 + 10) z' 1 

60 700 

+ - 2~ 1/3 e~^ /3 a s (281 a 3 + 10440) z~ 5/3 

1134000 V ; 

' 2 1/3 e i2w/3 a 2 (73769 a 3 + 6624900) z" 7/3 , \z\ > 1, (17) 



2619540000 



to be compared with those given in Refs. [3lllfl9ll2Q|[27] . According to Eq. (j2J) 
the z-zeros of i^W(z) are given approximately by 

Zs ~ e - in/2 (u + 9 s u 1/3 ), H»l: s = l,2,3,.... (18 



The usefulness of Eqs. ( ITT]) and (fTBl is restricted to the first values of s = 
1, 2, 3, . . ., due to the fact that \a s \ increases with s. This makes the the omit- 
ted terms in the expansion ( fTTl) to increase and invalidates their suppression. 
Besides this, \8 S \, given by (ED , increases with s and, for a given u, may reach a 
value such that expansion ffl3|) ceases from being useful. If zeros corresponding 
to a large s are to be considered, approximate values can be obtained follow- 
ing a procedure, already used by Magnus and Kotin [25] and by Cochran [3], 
consisting in taking logarithms in Eq. PJ), 



2 i/ log(z/2) + 1 2nvr = \og(T(u)) - log(-r(-i/)) 

+ 1 o4£i^^)- lo 4£i^7TX^V n = 0,±l,±2,..., (19) 

and approximating the logarithms of the gamma functions by their asymptotic 
expansions [U Eq. 6.1.40]. One obtains in this way 
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log(z/2) ~ log v — 1 — in/2 — i(n — j)tt/u + E 



OO 75 



2i» 



,a=i 2m(2m - l)u 



u 2m 

21/ v H^a-^J H^a+^k 



(20) 



Assuming |z| <C |z/| and keeping the dominant terms, it comes out for the 
z-zeros of K v (z) 



z n ~ e" irr 



12 2 i/ exp (-1 - i(n - . (21) 



Obviously, since we are considering values of z/ in the first quadrant, only large 
positive values of the integer n are compatible with the assumption \z\ <C \u\. 
Moreover, for a given value of \u\, the approximation (|2T|) looses accuracy 
as argz/ goes from it/ 2 towards and becomes completely useless for real v. 
A numerical comparison of the values obtained with Eqs. (fl8l) and ([21]) for 
moderate values of s and n shows that both labels coincide. 



4 Zeros of K v {z) for intermediate values of \v 



We have used a numerical procedure to explore the behaviour of the z-zeros 
of K u (z), as a function of u, for moderate values of \v\. In order to simplify 
the presentation of the results, two different situations have been considered: 
(i) variable \u\ with fixed argz/, and (ii) variable argz/ keeping \u\ fixed. The 
resulting trajectories of the first, second, and third zeros are represented in 
Figures 1, 2, and 3, respectively. For the evaluation of K v (z), expression 
has been used. More sophisticated algorithms can be found in the literature 
[2ll9lll3lll4f 1511161117(2 9 30 31 32j. but either they have been conceived only for 
particular (real or pure imaginary) values of the variable or the order, or 
they become unnecessarily complicated for our purpose of giving a merely 
qualitative description of the zeros. For the location of these, the Newton 
method has been used. Let us mention, however, other existing techniques 
like, for instance, the auxiliary tables for the evaluation of the z/-zeros of 
K u (z), for z real or pure imaginary, proposed by Cochran and Hoffspiegel [I], 
the finite element approximation method, applied by Leung and Ghaderpanah 
[22] to a very precise determination of the zeros of K n (z), and a procedure, 
applied by Segura [28] to the computation of zeros of Bessel and other special 
functions, that uses fixed point iterations and does not require the evaluation 
of the functions. 
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Im z 




Fig. 1. Trajectories of the first zero, z±, of K u (z). The continuous lines represent 
the trajectories followed by z\ as v varies whereas argf remains fixed and equal 
to m.7r/20, with m = (first trajectory from the left), 1, 2, . . . , 10 (horizontal tra- 
jectory). The continuous lines turn into short-dashed ones as the zero leaves the 
Riemann sheet — ir < arg z < ir. The trajectory corresponding to real v (m = 0) 
goes to (+oo, vr/4), with argzi — > —2ir, for v = 1/2. The nearly circular dashed arcs 
represent the trajectories of z\ for fixed integer \v\ and varying argf going from 
to 7r/2. The trajectories corresponding to \v\ = 5 and 10 are explicitly indicated. 

A double precision FORTRAN code was used to compute the expression 



£u(z) 



™ {z 2 /4) k 



(z/2) 



2i> 



ra- 



E 



> 2 /4) ^ 



(22) 



whose zeros coincide with those of the right hand side of Eq. ([3]). The two series 
were summed term by term until the absolute value of the ratio of the last 
term to the partial sum turned out to be less than 10 -18 . A warning message 
was foreseen for the case of the number of summed terms becoming larger than 
100, but this limit was never reached. Fortunately, zeros at large values of \z\ 
occur for values of v either of large modulus or in the neighbourhood of n + |, 
n — 0, 1, 2, .... In the first case, the series in (|22|) converge reasonably, from the 
numerical point of view; in the second case, analytical methods can be used to 
approximate the zeros. Typical runs of our procedure can be seen in Table 1, 
which shows intermediate results in the determination of the first, second and 
third zeros of K v (z), with \v\ = 21, aigv = 7tc/20. Initial approximate values 
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Fig. 2. Trajectories of the second zero, z%, of K u {z). The same convention as in Fig. 
1 has been adopted to represent the behaviour of z-2- The trajectory for real v goes 
to (+oo,37r/4), with arg^2 — > — 2tt, for v = 3/2. 



of the zeros, for a given v, were obtained by extrapolation of the trajectories 
followed by the zeros as v varies. Besides the smallness of the values of the 
real and imaginary parts of E v {z), their changes of sign in the successive steps 
of the Newton method make us to be confident about our results. 

Our figures show that for large \v\, with argz/ G (0, 7r/2), the zeros lie far 
from the origin in the half plane — it < axgz < 0, according to Eq. (|T8|) . 
When \v\ decreases, they approach the origin but, before reaching it, they 
cross the semiaxis argz = — it, leaving in this way the principal Riemann 
sheet —Ti < argz < n, for a certain |i/ cr it, s | whose value depends on argi/ and 
on the order s of the zero z s . The smaller aigu (between and vr/2) and the 
higher the label s of the zero, the larger the value of |f cr it,s|- ( in the case of 
real v, it can be shown that ^ cr i tiS = 2s — 1/2.) Therefore, except for the case 
of pure imaginary u, the number of zeros of K v {z) in the principal Riemann 
sheet is finite: only those zeros with label s — 1, 2, . . . such that |f cr it,s| < IH> 
for the corresponding value of argz/, remain in that principal Riemann sheet. 
If \u\ decreases again, being argz/ ^ 0, the zeros approach the origin spiraling 
clockwise infinitely, according to Eqs. (fTTj) and (|T2|) . The case of real v deserves 
a particular consideration. 
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-10 

-10 -5 5 Re 1 10 

Fig. 3. Trajectories of the third zero, 2:3, of K u (z). The lines represent the behaviour 
of Z3 with the same convention as in Fig. 1. The trajectory for real v goes to 
(+oo,57r/4), with arg^3 — > — 2ir, for v = 5/2. 

5 Zeros of K„(z) for real v. 

Let us now concentrate on the particular case of v being real. According to 
the second of Eqs. ([T]), the zeros of K v {z) appear in complex conjugate pairs. 
We refer, in what follows, only to those zeros of negative argument. Let us also 
recall that, for half-integer u, v — n + |, the Macdonald function reduces to 
a factor, different from zero in the finite plane, times a polynomial of degree 
n. Therefore, K v {z) presents exactly n z-zeros in each one of the Riemann 
sheets. Of course, the locations of these zeros are the same in all sheets. 

The aproximate expression (JT8l) for the first zeros of K v (z) [y ^> 1) is valid also 
in this case. Equation (ED), instead, is not applicable because the assumption 
\z\ \v\, used in its derivation, is not justified. Also, Eqs. (fTTll and (fT2ll are 
no more valid due to the possibility of \z\ becoming infinite for certain values 
of z/, a fact that invalidates Eq. ([5]), from which Eqs. ffTTT) and f|T2|) stem. 

According to Eq. (2), the z- zeros of K u (z) are obtained from the z-zeros of 
H^'(z) by a rotation of —tt/2 in the complex 2-plane. It is then immediate to 
deduce the behaviour of the zeros of K v (z) for varying real v from a previous 
discussion [5|6] of the zeros of H^\z). For convenience of the reader, we recall 
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Table 1 

Typical output of the Newton method applied to the location of the first, second 
and third zeros of K u (z), with \v\ = 21, aigv = 7n/2Q. Besides the successive 
approximations to the zeros, the values of S v (z), given by Eq. (|22p . are shown. 



Ski 


3zi 






14.00000000 


- 8.60000000 


- .127330E+01 


- .817656E-01 


14.02983390 


- 8.68222831 


.154771E+00 


.778107E-01 


14.02406924 


- 8.67520518 


.105223E-01 


.979004E-03 


14.02389450 


- 8.67463545 


- .594713E-04 


.105606E-04 


14.02389461 


- 8.67463884 


.429746E-07 


- .440585E-08 










10.90000000 


- 8.00000000 


- .151232E+01 


- .835798E+00 


10.98486986 


- 7.94986297 


- .232404E+00 


.536993E-01 


11.00073780 


- 7.95926858 


.207686E-01 


- .296534E-01 


10.99984823 


- 7.95696547 


.191900E-03 


- .214639E-03 


10.99983888 


- 7.95694790 


- .278481E-06 


.649298E-06 


10.99983889 


- 7.95694795 


.698624E-11 


- .156303E-10 






R£Jz 3 ) 




8.50000000 


- 7.00000000 


.741241E+00 


.345969E+01 


8.60000000 


- 6.92107409 


.978987E+00 


.293843E+01 


8.70000000 


- 7.02107409 


.658986E+00 


.247911E+01 


8.80000000 


— 7.12107409 


.598150E+00 


.174818E+01 


8.88169635 


- 7.22107409 


.831925E+00 


.871598E+00 


8.80046094 


- 7.31734050 


- .282066E+00 


.243279E+00 


8.82858965 


- 7.32026067 


.229670E-01 


.745812E-01 


8.82891727 


- 7.32660986 


.257068E-04 


- .695341E-03 


8.82889674 


- 7.32655881 


- .673649E-06 


- .726223E-05 


8.82889659 


- 7.32655825 


.622240E-10 


.575485E-09 



here, without demonstration, the main results of that discussion translated to 
those zeros in the lower half plane < 0. To be specific, let us focus on the 
trajectory followed by a zero, z s , as v decreases from +oo to 0. 

For very large v, z s is given approximately by Eq. (fT8|) . It lies far from the 
origin in the quadrant — n < argz < —n/2. As v decreases, it moves upwards 
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and to the right, reaching the horizontal semiaxis aigz = — ix for v = 2s — 
1/2. Then, it makes a nearly circular quarter of a turn, crosses the vertical 
semiaxis argz = — 37r/2 for v = s — 1/3, and goes nearly horizontally towards 
(+oo, (2s — l)ir/4) as v — > s — 1/2 from above. 

The behaviour of z s for v < s — 1/2 is not uniquely defined. In fact, the implicit 
function z s {v) defined by the condition K v {z) = presents a logarithmic 
branch point at v = s — 1/2, and the values of z s for v < s — 1/2 depend 
on the chosen branch. One can pass continuously from the values of z s for 
v > s — 1/2 to those for v < s — 1/2 if one avoids the branch point by the 
usual procedure of adding a small imaginary part to v as its real part crosses 
the value s — 1/2. But the resulting values of z s for v < s — 1/2 depend on 
the sign of that imaginary part. Since this ambiguity in the values of z s has 
been thoroughly discussed in Ref. [6], in the context of the z-zeros of H^(z), 
we do not consider necessary to pursue further. 

Besides those described above, K„(z) presents an infinite set of zeros near 
and at the left of the vertical semiaxis aigz = — 37r/2 for every integer v = 
n [2H p. 513]. As v increases, they approach that semiaxis, cross it when 
v = n + 1/3, and move nearly horizontally to the right, going to infinity as 
v — » n + 1/2, which is a logarithmic branch point for the position of the zeros 
as a function of u, as it was shown in Refs. [5|6] for the zeros of HjP(z). 
The zeros make a discontinuous jump of ±z7r/2, according to how the branch 
point is circumvented, move nearly horizontally to the left, reach the vertical 
semiaxis argz = — 37r/2 when v = n + 2/3, and go to the left, to end, for 
v = n + 1, at positions intermediate between those that they occupied for 
v = n. A figure illustrating that behaviour can by obtained by rotating Fig. 1 
in ref. [5] by an angle of — 7r/2. 
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